On vacuum gravitational collapse in nine dimensions 
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We consider the vacuum gravitational collapse for cohomogeneity-two solutions of the nine dimen- 
sional Einstein equations. Using combined numerical and analytical methods we give evidence that 
within this model the Schwarzschild-Tangherlini black hole is asymptotically stable. In addition, we 
briefly discuss the critical behavior at the threshold of black hole formation. 



Introduction and setup: Over the past few years there has been a surge of interest in higher dimensional gravity, 
motivated by several reasons. The main reason comes from the new brane-world scenario in string theory according 
to which we live on a three dimensional surface (called a brane) in a higher dimensional spacetime In order to 
understand the phenomenology of this scenario and its experimental predictions, like production of black holes in 
the next generation of colliders, it becomes important to study solutions of Einstein's equations in more than four 
dimensions. Another reason, closer to our own motivation, has nothing to do with string theory and comes from general 
relativity itself. Viewing the dimension of a spacetime as a parameter of the theory is helpful in understanding which 
features of general relativity depend crucially on our world being four dimensional and which ones hold in general. 
Last, but not least, extra dimensions are fun. 

In a recent paper some of us showed that in five spacetime dimensions one can perform a consistent 
cohomogeneity-two symmetry reduction of the vacuum Einstein equations which ~ in contrast to the spherically 
symmetric reduction ~ admits time dependent asymptotically flat solutions. The key idea was to modify the standard 
spherically symmetric ansatz by replacing the round metric on the three-sphere with the homogeneously squashed 
metric, thereby breaking the SO{A) isometry to 5*0(3) x U{1). In this way the squashing parameter becomes a 
dynamical degree of freedom and the BirkhofF theorem is evaded. This model provides a simple theoretical setting for 
studying the dynamics of gravitational collapse in vacuum, both numerically 2] and analytically ■ 

As mentioned in , similar models can be formulated in higher D = n + 2 dimensions as long as the corresponding 
sphere 5" admits a non-round homogeneous metric, i. e. there exists a proper subgroup of the orthogonal group 
SO{n + 1) which acts transitively on 5". According to the classification given by Besse such transitive actions 
exist on all odd dimensional spheres. For example, the group SU{n -\- 1) acts transitively on S*^""*"^ and the group 
Sp{n -\- 1) acts transitively on 5"*"+'^. It is natural to ask whether the properties of gravitational collapse found in [2| 
are typical for this class of models or whether new phenomena appear in higher dimensions. A systematic analysis of 
this question appears hopeless in view of the fact that the number of degrees of freedom (squashing parameters) grows 
quickly with dimension, thus it seems useful to look at specific examples to get better understanding. In this note 
we consider a model with one degree of freedom which describes the squashing of the seven sphere. More concretely, 
we regard S"^ as the coset manifold Sp{2) / Sp{l) ~ SO{5)/ SO{3), or equivalently, as the bundle over the S'^ base 
space with 50(5) x 50(3) invariant metric. In the past, such the squashed seven sphere has attracted a great deal 
of attention in the context of eleven dimensional supergravity 

We parametrize the metric in the following way 

ds^ = -Ae'^^dt^ -f A-^dr"^ + r^dft^, (1) 
where dfl^ is the metric on the unit squashed 5"^ 6| 

dnj = le'^id^.' + e^) + le-'^Ef, (2) 



B, A, and S are functions of t and 



i sin ^ (ffi - CTi), E, = cos^ ^ + ^"^^ ^ ^i' (3) 



and fji and ai are two sets of left-invariant one-forms on SU (2) 

CTi = cos 4'dO -I- sin ■0 sin 9d(t), a2 = — sin tpdO -\- cos tp sin 6d(j), 1T3 = di/j -I- cos 6d(j), (4) 

with ai given by identical expressions in terms of {ip,9,(j)). In this ansatz the S'0(8) isometry of the round 5^ is 
broken to 50(5) x SU{2). This is a special case of a more general ansatz for which the fibers are allowed to be 
squashed themselves P]. 



2 



Substituting the ansatz into the vacuum Einstein equations, Rap = 0, we get equations of motion for the 
functions A{t,r),5{t^r) and B{t,r) (in the foUowing we use overdots and primes to denote dt and dr, respectively) 

A' = + 1 (48e-3B _ ^2e-ios + Ge^^) - 3r (^e^'A-^B^ + AB'') (5) 
A = -6rABB', (6) 
6' = -3r (^e^^A-^B^ + S'^) , (7) 

(e*A- V^b) ■ - {e-'Ar'B')' + ^e^V^ (Ge^^s - e^^ - Se-^"^) = 0. (8) 

It foUows immediately from the above equations that if S = (no squashing), then 6 — and either A = 1 
(Minkowski), or A = 1 — rf^/r^ (Schwarzschild-Tangherlini), where rh is the radius of the black hole horizon. We 
showed in 0] that the analogous two solutions in five dimensions play the role of generic attractors in the evolution 
of regular asymptotically flat initial data (small and large ones, respectively) and the transition between these two 
outcomes of evolution exhibits the type II discretely self-similar critical behavior. The aim of this note is to show 
that these properties are also present in nine dimensions which suggests that they are general features of vacuum 
gravitational collapse for this class of models. 

Linear stability and quasinormal modes: Before presenting numerical evidence for the nonlinear stability of the 
Schwarzschild-Tangherlini solution, we want to discuss the results of linear perturbation theory. Linearizing equations 
©-(IHl) around the Schwarzschild-Tangherlini solution we obtain the linear wave equation for the perturbation SB{t, r) 

6B ~ ^Aoir'AoSB'Y + ^^SB = 0, Ao = l-^, (9) 

where we have used the scaling freedom to set the radius of the horizon rh — 1- Introducing the tortoise coordinate 
X defined by dx/dr — Aq ^, and substituting 5B{t,r) = e~*'''*r~''/^it(a:) into (PJ we get the Schrodinger equation on 
the real line — cx3 < a; < oo 

d'^'^ T.w ^^ ,2 T.. ^ 1^ l\/99 49\ 

-^+yi^i.^))^-k^^^ V^^) = l[^--,)[y + ^)- (10) 

This equation is a special case (corresponding to the gravitational tensor perturbation with Z = 2) of the master 
equation for general perturbations of the higher dimensional Schwarzschild-Tangherlini solutions derived independently 
by Gibbons and HartnoU @ and Ishibashi and Kodama Since the potential in (|10|l is everywhere positive, there 
are no bound states, which implies that the Schwarzschild-Tangherlini black hole is linearly stable. 

We computed quasinormal modes, i.e. solutions of equation H10() sat isfy in g the outgoing wave boundary conditions 
u e^**^^ for X icx), using Leaver's method of continued fractions [iflllll- Substituting 

= ^^"■E«»(^) ' = (11) 

into pO(l we get a 12-term recurrence relation for the coefficients a„. The prefactor in fulfills the outgoing wave 
boundary conditions both at the horizon and at infinity, hence the quasinormal frequencies are given by the discrete 
values of k for which the series in converges at r = c». Using Gaussian elimination we reduced the 12-term 
recurrence relation to a 3-term one. According to Leaver [lfil |. the convergence of the series X]^o equivalent 
to a condition under which the 3-term recurrence relation has a minimal solution. Such the condition, formulated in 
terms of continued fractions, yields the transcendental eigenvalue equation which has to be solved numerically. The 
full details of this computation will be described elsewhere. Here, in order to interpret numerical results we need only 
to know that the fundamental quasinormal mode has the frequency = 3.4488 — 0.8601i (in units t~^^^- This mode 
is expected to dominate the process of ringdown. The higher modes are damped much faster (for example, the first 
overtone has the frequency k\ ~ 2.7548 — 2.6116i) so in practice they play no role in the dynamics. 

Numerical results: Using the same finite- difference code as in we solved equations (5)-(8) numerically for 
several families of regular initial data. We found that the overall picture of dynamics of gravitational collapse looks 
qualitatively the same as in five dimensions, that is, we have dispersion to the Minkowski spacetime for small data and 
collapse to the Schwarzschild-Tangherlini black hole for large data (see Fig. la). In the latter case we looked in more 
detail at the asymptotics of this process. We found that at some intermediate times the solution settling down to 
the Schwarzschild-Tangherlini black hole is well approximated outside the horizon by the least damped quasinormal 
mode. This is shown in Fig. lb. 
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(a) Formation of a black hole 



(b) Quasinormal ringing 



FIG. 1: (a) For a black-hole solution we plot a series of late-time snapshots of the mass function m{t,r). The initial mass 
function (solid line) has the form of a kink asymptoting the total mass m(oo) = 0.575. During the evolution the horizon 
develops at = 0.606 which corresponds to the black hole mass msH = rf^ = 0.05. Thus, only a small fraction of the total 
mass gets trapped inside the horizon while the remaining mass is being radiated away to spatial infinity, as is clearly seen from 
the plot, (b) We plot the time series In \B(t,ro)\ at ro = 5 for the same solution as in (a), and superimpose the fundamental 
quasinormal mode with frequency fco = (3.4488 — 0.8601 i)/r(,. On the time interval 10 < i < 16 we get very good agreement 
between these two curves which confirms a well-known fact that quasinormal modes encode an intermediate time behavior of 
solutions at a fixed point in space [l^. 
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(b) Black-hole mass scaling 



FIG. 2: (a) For a near-critical evolution we plot the function B at some late central proper time To and superimpose the next 
three echoes (at times r„) shifted by ln(r) — > ln(r) +nA. Minimization of the discrepancy between the profiles yields A ~ 0.78. 
(b) For supercritical solutions the logarithm of black hole mass tubh is plotted versus the logarithmic distance to criticality. 
A fit of the power law ln(m_B£r) ~ "fln{p — p") + const yields 7 ~ 1.64. The small wiggles around the linear fit are the imprints 
of discrete self-similarity; their period is equal to 2.87 in agreement with the theoretical prediction 6A/7. 



At the threshold for black hole formation we observe the type II discretely self-similar critical behavior with the 
echoing period A « 0.78 and the black-hole mass scaling law with the universal exponent 7 « 1.64. This is shown 
in Fig. 2. Note that in eight space dimensions, mass has the dimension of length^, hence 7 = 6/A, where A is the 
eigenvalue of the growing mode of the critical solution. Strangely enough, the product AA « 2.85 is approximately 
the same (up to numerical errors) as in five dimensions 0. We do not know whether there is any deeper meaning 
behind this numerical coincidence. 
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Final remarks: In this note we have focused on similarities between the models in five and nine dimensions, however 
at the end we would like to mention two interesting qualitative differences which in our opinion are worth investigation. 

The first difference is the existence of the second static solution corresponding to the non-round homogeneous 
Einstein metric on S'^ 

B = ^, A = 9-5''"/'[l-'iy 5 = 0. (12) 

This solution is asymptotically conical so it does not participate in the dynamics of asymptotically flat initial data. 
It would be interesting to study the evolution of non-asymptotically flat initial data and determine a dynamical role 
of the solution (|12|l . as well as other known explicit solutions, like for example the Spin(7) solution [l3l| . 

The second difference is the lack of monotonicity of the mass function m{t, r), defined by ^ = 1 — m{t, r)/r^ . From 
the hamiltonian constraint (5) we obtain 

m' = 3/ (e^'A-^B^ + AS'') + ^r^ (42 - 48e-3B + la^-ios _ q^^b^ . (13) 

The potential term on the right hand side of this equation has the local minimum equal to zero at i? = 0, the local 
maximum at B = In 5/7, and is negative for large positive values of B. Thus, the mass density may be locally negative 
and, indeed, it is easy to construct initial data with large regions of negative mass density. This suggests a possibility 
of violating the weak cosmic censorship (note that monotonicity of mass is essential in the Dafermos and Holzegel 
proof 0). Although our preliminary numerical attempts failed to produce a counterexample, i.e. a generic naked 
singularity, this problem deserves more systematic investigation. We should stress that despite locally negative mass 
density the total mass m(oo) is guaranteed to be nonnegative by the positive mass theorem in higher dimensions |14| . 
It would be interesting to prove this fact in our model in an elementary manner. In particular, for time-symmetric 
initial data it follows from equation (|13|l and the requirement of regularity at the origin, m(0) = 0, that 

m(oo) = ^'^^^^""""^ efo^P'"'"'"' (^r'^B'^ -f ir'^(42 - A^e'^^ + 12e-i°^ - Ge^^)) dr. (14) 

We challenge the readers to show that this integral is nonnegative for all "reasonable" functions B{r) which vanish 
at the origin and have compact support (or fall off faster than at infinity). Finally, we remark that there is 
an analogous problem for the black-hole boundary condition, m{rh) = where one wants to prove the Penrose 
inequality m{oo) > m{rh). 
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